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ABSTRACT 


This  report  describes  a  parametric  study  of  stiffened  cones  and  cylin¬ 
ders  to  analytically  determine  the  effects  of  a  longitudinal  steady 
state  driving  force  on  the  radial  velocity  of  axially  symmetric  shells. 
The  "lumped-mass"  technique  is  utilized  to  predict  the  natural  fre¬ 
quencies,  mode  shapes  and  impedances  (point  and  transfer)  for  eight 
cases.  Although  the  natural  frequencies  are  not  too  sensitive  to 
changes  in  shell  thickness  and  frame  area,  the  transfer  impedances 
(hence  radial  velocities)  vary  significantly  for  the  cases  studied. 
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NOMENCLATURE 
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M- 
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X1,  f,  Zl 

■t 

as 

CP 
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lumped  mass  at  joint  j 
deflection  at  joint  j 
stiffness  of  each  member  (p) 
area  of  member  (p) 

Young's  Modulus 

Poisson's  ratio 

length  of  member  (p) 

moments  of  inertia  of  member  (p) 

polar  moment  of  inertia  of  member  (p) 

shear  distribution  factors 

Joint  (j)  restraint  factor 

coordinates  normal  to  longitudinal  axis 

longitudinal  coordinate 

skewed  coordinate  systems  for  member  (p) 

time  in  seconds 

eigenvalues 

eigenvectors 

natural  frequency  (f  =  ■—  cycles  per  second), 
generalized  mass 
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INTRODUCTION 


The  coupling  of  shell  modes  is  defined  as  the  interaction  of  the 
shell's  stretching  modes  and  the  shell's  lobar  modes  in  a  vibrating 
shell.  This  interaction  affects  the  sound  radiation  from  the  shells 
when  they  are  subjected  to  steady  state  driving  forces.  Thus,  a 
cylinder  being  driven  along  its  axis  not  only  creates  sound  pressure 
in  the  direction  of  the  force  but  also  in  the  'radial  direction.  In 
this  report,  the  mechanical  coupling  of  the  modes  is  investigated 
but  not  the  resulting  sound  radiation. 

Lobar  frequencies  of  shells  are  usually  calculated  separately  from 
the  stretching  frequencies  for  the  sake  of  simplifying  the  calcula¬ 
tions.  But  to  determine  the  true  interaction  between  the  two,  they 
must  be  considered  simultaneously.  To  facilitate  the  latter,  a 
computer  program  was  developed  [  1,  2,  3,  &  4]*  that  can  calculate 
the  frequencies,  mode  shapes,  and  mechanical  impedances  of  a  shell 
of  any  (arbitrary)  geometry.  The  technique  used  is  commonly  called 
the  "lumped-mass"  method.  By  this  is  meant  that  structural  systems, 
in  this  case  stiffened  cylinders  and  cones,  are  idealized  as  a  grid- 
work  of  plates  and  beams  whose  masses  are  concentrated  at  their 
respective  centers  of  gravities.  In  the  case  of  a  stiffened  cylinder, 
the  gridwork  consists  of  a  series  of  flat  plates  (Pig.  l)  related  by 
direction  cosines  to  form  the  shell,  and  straight  beam  members,  also 
related  by  direction  cosines,  to  form  the  stiffeners.  Continuity 
and  equilibrium  are  satisfied  at  the  joints  between  members.  Thus 
the  stiffness  (and  flexibility)  of  the  entire  structure  can  be 
determined,  and  having  established  the  mass  distribution,  the  equations 
of  motions  can  be  written  and  solved  to  determine  natural  frequencies 
and  mode  shapes.  The  point  and  transfer  impedances  can  then  be 
determined. 


*  Refers  to  references  (pg  18) 
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In  this  study,  the  objective  is*  to  determine  the  natural  frequencies, 
mode  shapes,  and  Impedances  of  two  basic  shells  of  revolution; 
cylinders  and  cones.  After  choosing  the  diameters  and  lengths, 
several  combinations  of  shell  thicknesses  and  stiffener  areas  were 
considered.  The  results  indicate  to  what  extent  frame  size  and  shell 
thickness  influence  the  vibrational  behavior  of  cylinders  and  cones, 
which  in  turn  influences  the  nature  and  strength  of  the  sound  radia¬ 
tion  from  these  shells.  Ultimately,  a  thorough  knowledge  of  the 
effects  of  the  shell  and  stiffener  geometry  on  the  vibration  and 
radiation  characteristics  can  lead  to  structures  with  minimal  adverse 
behavior. 


6 


II 

THEORETICAL  DEVELOPMENT 

The  initial  operation  in  analyzing  stiffened  shells  (Pig.  la)' by  the 
lumped  mass  technique  is  to  idealize  the  shell  and  stiffeners  by  a 
system  of  flat  plates  and  bars  (Pig.  lb).  Hrennlkoff [5] has  shown 
that  the  flexibility  of  this  system  can  be  duplicated  by  using  a 
mathematically  equivalent  framework  of  elastic  bars  (Pig.  lc).  This 
equivalent  framework  is  then  used  as  the  mathematical  model  for  the 
structure . 

Each  member  of  the  framework  has  typical  beam  properties,  but  these 
properties  are  mathematical  equivalents  only,  and  do  not  have  direct 
physical  meaning.  A  typical  cross  section  and  plan  view  of  an 
unstiffened  shell  are  shown  in  Pigs.  2a  and  2b. 

The  properties  of  the  members  in  these  panels  are: 


The  addition  of  a  stiffener  (Pigs.  3a  and  3b)  contributes  to  the 
stiffness  of  the  circumferential  frames. 

The  only  properties  which  change  are  and  I 


Als  A°  +  A( frame) 
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^xls  Ixl  +  Ix( frame) 

Iyl  =  Iy( frame) 

Since  the  amount  of  damping  inherent  to  thin  shells  is  very  small 
compared  with  the  value  for  critical  damping,  the  effects  on  the 
natural  frequencies  and  mode  shapes  are  negligible.  Thus,  the  only 
parameters  needed  to  determine  the  natural  frequencies  are  the 
structural  geometry,  mass  distribution  and  end  conditions,  which  are 
all  determined  in  establishing  the  equivalent  gridwork  of  elastic 
beams . 


The  equations  of  motion  for  free  vibrations  are: 

• '  n 

m.  D.  +2  S..  D.  =  0  (l) 

J  J  i=l  1 

where  j  is  the  joint  number. 

Each  joint  (j)  is  common  to  two  or  more  members  (p),  and  the  joint 

stiffness  (Sjj)  is  comprised  of  the  sum  of  the  stiffnesses  of  each 

member  (S  ).  The  stiffness  of  each  member  at  the  centroid  of  the 
P 

member  is: 


S 


P 


-0  1 


L  o- 


ap.  a 


(1.1) 


where 


aa8  =12  EIxA/[i3A  +  24(1  +  M-  ] 

a33  =12  EIyA/  [i3A  +  24(1  +U  )KyiIy)  ] 
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a44  -  JG/(1  -  t)J 
asB  -  Bl/i 
a*8  =  Eiyi 

Before  the  stiffnesses  of  each  member  (S  )  are  summed  to  get  the 
joint  stiffness  (Sjj),  each  member's  stiffness  Is  rotated  and  trans¬ 
lated  from  Its  own  coordinate  system  to  the  coordinate  system  of  a 
common  origin. 


Thust 

m  rn 

Sp°  -  ITBS  BAL- 

where 

L  — 


Ki  2  0 

0  K,a 


' cosx1 X 
cosy1  x 
cosz1  x 


cosx1  y 
cosy1  y 
cosz1  y 


cosx1  z 
cosy1  z 
cosz1  z 


(1.2) 


(1-3) 


1 - 

0  1 

0  0  1 

o  (zP-z0)  (yc-yp)  1 

(zfj-zp)  0  (xp-Xo)  o 

(yp-y0)  (xo-xp)  o  o 


o 


1  j 

0  1 


(1.4) 


Once  the  stiffnesses  of  all  the  joints  are  calculated  at  the  origin, 
they  can  be  used  to  form  the  stiffness  matrix  and  the  ^2-^ ^ j  1  ^1 
matrix  In  the  equation  of  motion. 


n 

2 

1-1 


n 

2 

k-1 


K  /  j 


(2) 
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Thus,  in(Fig.  4  ) 


n 


2  S 
i=l 


o  W3 

ji  1 


[Si°+  si+  s°] 


-s 


o 

1 


-S? 

s?+  s£  +  s£ 
+s°+  s°e 


0 

-S| 


etc  .■ 


Or,  in  matrix  form 


2SOD°«£PD0  (2.2) 

i=l 

where  S°  is  the  n  x  n  stiffness  matrix  at  the  origin 
and  I f  is  the  n  x  1  deflection  vector  at  the  origin 


Before  the  natural  frequencies  of  the  system  are  calculated,  the 
stiffness  matrix  at  the  origin  (S°)  is  inverted  to  form  an  influence 
coefficient  matrix  (6°). 

[s0]"1  =  5°  (2.3) 

The  influence  matrix  forms  the  characteristic  matrix  and  is  transferred 
back  to  the  initial  coordinate  system  at  the  Joints. 

6  >=  BTLTFT6°LBF  (2.4) 

F  is  a  force  matrix  with  unit  forces  at  each  mass  for  each  degree  of 
freedom. 

The  equations  of  motion  then  can  be  written  in  terms  of  the  influence 
coefficient  matrix: 

D  -  -  6  M  D  (3) 
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where : 


D  is  the  displacement  vector  (6n  x  l) 

D  is  the  acceleration  vector  (6n  x  l) 

6  is  the  flexibility  matrix  (6n  x  6n) 

M  is  the  diagonal  mass  matrix  ( 6n  x  6n) 

The  natural  frequencies  can  be  determined  by  assuming  a  peroldic 
displacement 

D  =  qpsinot  (3.1) 

and  inserting  into  equation  (3)  to  get  the  frequency  equations 

i  -  6M  «  0  (4) 

J 

A  more  expedient  form  of  equation  (4)  is  obtained  if  the  mass  matrix 
is  factored  so  that: 

M  =  ATA  (4.1) 


where : 


Premultiplying  equation  (4)  by  A  and  using  (4.1)  results  in: 


I  T 

-  -  A6A  =  0 

of 


(5) 
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Since  X6 X  is  a  symmetrical  matrix,  computer  storage  for  only  half 
of  the  non-diagonal  terms  is  required. 

The  solution  of  equation  (5)  is  accomplished  by  using  the  Modified 
Givens  Method  jg]  to  determine  the  eigenvalues  (o>)  and  eigen- 

vectors  (9).  Due  to  the  operation  whereby  X6X  is  used,  the  eigen-  ^ 

vectors  are  not  the  true  ones,  but  are  reorientated.  The  true  eigen¬ 
vectors  are  then  obtained  by:  t 

<p  *  X-1  8  (6) 

The  equation  for  the  undamped  transfer  and  point  impedances,  using 
the  eigenvalues  (to)  and  eigenvectors  (cp)  of  the  stiffened  shell,  is 


ZJk  = 


&  n 

M  2 

r=l  ofn 


CD 


, where  j  =  k  for  point  impedance. 


(7) 


i 
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Ill 

DESCRIPTION  OF  MATHEMATICAL  MODELS 

Eight  models  were  analyzed,  four  cones  and  four  cylinders  (Pig.  5 ) • 

For  all  cases,  the  end  conditions  were  the  same;  zero  radial  deflection, 
zero  longitudinal  slope  and  free  longitudinally  at  the  end  with  the 
steady  state  driving  force  (Fe^ffi^).  The  length  and  diameter  of  the 
cylinder  were  held  constant,  and  the  shell  thickness  and  area  of 
frames  were  varied.  Similarly,  in  the  cone  the  length  and  diameter 
were  held  constant  and  the  shell  thickness  and  areas  of  frames  were 
varied. 

Since  four  cases  were  studied  for  each  the  cone  and  cylinder,  six 
comparisons  can  be  made  for  each  geometry  (Table  l) .  In  both  the 
cone  and  cylinder,  two  shell  thicknesses  and  two  areas  of  frame  were 
considered.  Therefore,  the  cases  to  be  compared  are: 


3.1 


Comparisons 


TABLE  I 


a.  Doubling  the  frame  area  (l  In2  to  2  in3 )  while  keeping  shell 
thickness  constant  (1/4  in  ). 

b.  Doubling  the  frame  area  (l  in3  to  2  in3)  while  keeping  shell 
thickness  constant  (l/2  in). 

c.  Doubling  the  shell  thickness  ( 1/4  in  to  1/2  in)  while  keeping 
the  frame  area  constant  ( 1  in3 ) . 


d.  Doubling  the  shell  thickness  (1/4  in  to  1/2  in)  while  keeping 
the  frame  area  constant  ( 2  In3 ) . 

e.  Doubling  both  the  shell  thickness  ( 1/4  in  to  1/2  in)  and 
frame  area  ( 1  in3  to  2  in3  ) . 


f.  Doubling  the  shell  thickness  (1/4  in  to  1/2  in)  while  halving 
the  frame  area  ( 2  in3  to  1  in3  ) . 
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IV 

RESULTS 


«• 

The  most  significant  variations  In  results  for  the  cases  of  cones  and 
cylinders  were  in  the  impedances  (point  and  transfer) .  The  transfer  * 

impedances  are  the  most  important  to  this  study,  since  they  determine  j 

the  longitudinal  steady  state  driving  force  required  to  impart  a  unit 
velocity  to  either  the  shell  plating  or  frames  in  a  radial  direction. 

For  each  of  the  four  cases  of  cylinders  and  cones,  the  point  and 
transfer  impedances  were  determined  at  the  center  stiffener  and  at 
a  point  on  the  shell  midway  between  the  center  stiffener  and  the 
adjacent  stiffener  (Figs.  10  through  17)-  Si^ce  the  lower  lobar 
modes  (m  «  1,  n  *>  2,  3,  4)  are  the  most  efficient  sound  radiators, 
the  impedances  for  these  modes  are  most  significant.  Tables  4  and  5 
give  the  stiffnesses  (K  «  ia>Z)  and  Impedances  (point  and  transfer) 
for  the  lower  modes  for  the  various  cases  of  cylinders  and  cones. 

Table  6  compares  ratios  of  transfer  impedances  and  shows  that  the 

cones  and  cylinders  react  differently  to  the  parameter  changes.  I 

The  highest  increase  in  the  cylinders  is  derived  by  increasing 

the  shell  thickness  (comparison  c  &  d)  whereas  the  highest  increase 

in  the  cone  is  derived  by  increasing  the  frame  area  (comparison 

a  &  b)  . 

Contrary  to  the  impedance  results,  the  percent  changes  in  frequencies 
(stretching  and  lobar)  are  not  too  sensitive  to  the  changes  in  shell 
thickness  and  frame  areas.  For  each  case,  there  are  approximately 
25  frequencies  below  1  kc(l03cps).  The  frequencies  for  the  first 
four  longitudinal  modes  (m  «  1,  2,  3,  4)  and  first  six  lobar  modes 
(n  «■  2,  3,  4,  5,  6)  are  shown  in  Tables  2  and  3.  The  percent  changes 
in  frequencies  are  shown  in  Table  7. 
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For  both  the  cylinders  and  cones,  the  largest  increase  in  frequencies 
for  the  lowest  lobar  modes  (m  ■«  1,  n  •*  1,  2)  can  be  derived  by  in¬ 
creasing  the  shell  thickness  and  decreasing  the  area  of  frames 
(comparison  f  of  Table  1). 

The  analyses  of  the  cones  also  revealed  that  the  lobar  mode  shapes  in 
some  cases  are  "mixed"  (see  figures  6,  7 >  8,  9)*  In  some  cases  the 
frames  and  shells  have  different  mode  shapes  while  in  other  cases  the 
shell  and  frames  at  one  end  of  the  cone  have  different  mode  shapes 
from  that  of  the  shell  and  frames  at  the  other  end. 

The  percent  changes  in  frequency  for  the  lobar  modes  (n)  decrease 
with  longitudinal  modes  (m)  for  both  the  cones  and  cylinders.  The 
fourth  longitudinal  mode  (m  =  4)  has  the  smallest  change.  Since 
the  half-wave -length  for  this  mode  approximately  equals  the  length 
between  frames,  the  shell  behaves  somewhat  as  an  unstiffened  shell  and 
thus  the  changes  in  frequency  would  expectedly  be  small. 


V 


DISCUSSION  AND  CONCLUSION 

For  the  cases  studied,  the  variations  in  shell  thickness  and  frame  area 
have  their  greatest  effects  on  the  transfer  impedances.  Thus,  it  is. 
apparent  that  the  radial  velocity  of  both  shell  and  frames  can  be 
significantly  reduced  by  an  expedient  choice  of  frame  area  and  shell 
thickness.  The  effect  of  frame  spacing  was  not  investigated  in  this 
study,  but  it  would  also  have  a  significant  effect. 

Considering  the  two  parameters  studied  (shell  thickness  and  frame  area) 
the  most  efficient  way  of  raising  the  transfer  impedances  (thus  re¬ 
ducing  the  radial  velocity)  in  the  cylinder  is  to  keep  the  area  of 
frames  constant  while  increasing  the  shell  thickness.  In  the  cones, 
however,  the  opposite  is  true.  The  most  effective  means  of  increasing 
the  transfer  impedances  are  to  keep  the  shell  thickness  constant  while 
increasing  the  frame  area. 

Unlike  the  large  variation  in  impedances,  the  natural  frequencies 
(stretching  and  lobar)  do  not  substantially  change  for  all  the  cases 
of  cylinders  and  cones.  The  greatest  changes  are  accomplished  by 
increasing  the  shell  thickness  while  decreasing  the  frame  area. 

One  of  the  most  significant  results  is  the  fact  that  in  the  cones 
there  exist  "mixed"  lobar  modes.  These  "mixed"  modes  start  as  low 
as  the  fourth  mode  and  take  two  forms;  in  some  cases  the  mode  shapes 
at  either  end  of  the  cone  differ,  and  in  others  the  shell  has  one  mode 
shape  while  the  frames  have  another.  The  significance  of  these 
phenomena  is  that  only  an  analytic  method  that  does  not  predetermine 
a  mode  shape  (such  as  the  lumped  mass  methods)  will  determine  these  modes. 
Since  the  "lumped  mass"  technique  used  is  valid  for  any  geometry, 
symmetrical  as  well  as  not,  this  method  can  determine  these  "mixed" 
modes  with  equal  facility  as  well  as  the  regular  lobar  modes. 
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The  results  of  this  study  Indicate,  that  expedient  she.lJ  designs  cam 
significantly  reduce  the  radial  velocity  of  the  shells  when  subjected 
to  steady  state  longitudinal  driving  force. 
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(3B) 

FIGURE  3  CROSS  SECTION  AND  PLAN  VIEW  OF  STIFFENED  SHELL 
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POINT  IMPEDANCfc  -  RADIAL  STIFFENER  DRIVE 
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